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Abstract 

Verticillium wilt is a major threat to alfalfa (Medicago sativa) and many other crops. The model legume Medicago 
truncatula was used as a host for studying resistance and susceptibility to Verticillium albo-atrum. In addition to pre- 
senting well-established genetic resources, this wild plant species enables to investigate biodiversity of the response 
to the pathogen and putative crosstalk between disease and symbiosis. Symptom scoring after root inoculation and 
modelling of disease curves allowed assessing susceptibility levels in recombinant lines of three crosses between 
susceptible and resistant lines, in a core collection of 32 lines, and in mutants affected in symbiosis with rhizobia. 
A GFP-expressing V. albo-atrum strain was used to study colonization of susceptible plants. Symptoms and colo- 
nization pattern in infected M. truncatula plants were typical of Verticillium wilt. Three distinct major quantitative 
trait loci were identified using a multicross, multisite design, suggesting that simple genetic mechanisms appear to 
control Verticillium wilt resistance in M. truncatula lines A17 and DZA45.5. The disease functional parameters varied 
largely in lines of the core collection. This biodiversity with regard to disease response encourages the development 
of association genetics and ecological approaches. Several mutants of the resistant line, impaired in different steps 
of rhizobial symbiosis, were affected in their response to V. albo-atrum, which suggests that mechanisms involved in 
the establishment of symbiosis or disease might have some common regulatory control points. 

Key words: biodiversity, Medicago truncatula, nodulation mutants, partial resistance, quantitative trait loci, root disease, 
vascular wilt, Verticillium albo-atrum. 
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Introduction 

Vascular wilt diseases are major constraints for plant produc- 
tion. They are characterized by the colonization of the plant's 
xylem vessels by the pathogen and are mainly caused by soil- 
borne microorganisms, which might be bacteria or fungi, such 
as Fusarium or Verticillium (Agrios, 2005). The most important 
plant pathogenic Verticillium species are Verticillium dahliae 
and Verticillium albo-atrum, infecting more than 200 plant spe- 
cies and destroying high amounts of crops worldwide (Fradin 
and Thomma, 2006; Klosterman et al, 2009). The colonization 
of xylem vessels by spores and hyphae, together with gum for- 
mation by the plant cells, induce clogging of the vessels, stop 
water flow, and result in the typical wilting symptoms. At the 
end of the infection cycle, the fungus grows outside of the stele 
and colonizes other tissues (Agrios, 2005). Owing to survival 
structures, which remain viable for many years in the soil, and 
the protected localization in infected plants, Verticillium wilt is 
difficult to control. So far the most efficient way is by breed- 
ing resistant or tolerant varieties. Examples of polygenic and 
monogenic dominant resistance have been described for several 
plant species and resistance loci or genes have been identified. 
A well-studied example is the Ve locus in tomato which confers 
resistance against race 1 of V. dahliae and V albo-atrum (Fradin 
and Thomma, 2006). The tomato Ve locus has been cloned 
and shown to contain two genes encoding receptor-like pro- 
teins with extracellular leucine-rich repeat (LRR) domains, Vel 
and Ve2, which are each able to confer resistance to susceptible 
potato (Kawchuk et al, 2001). In Arabidopsis, the tomato Vel 
gene has been reported to confer resistance to V. dahliae and 
V. albo-atrum (Fradin et al, 2011). Studies on the interaction 
between lettuce and V. dahliae indicated that a single dominant 
gene with homology to tomato Vel and Ve2 controls resistance 
against race 1 of the fungus in this species (Hayes et al. , 201 1). 
Similarly, Ve homologues are proposed to be involved in toler- 
ance to V. dahliae in cotton and mint (Vining and Davis, 2009; 
Zhang et al , 201 1). The V dahliae effector Avel , recognized by 
tomato VE1, has been cloned recently (de Jonge et al., 2012). 

Because of the high protein content of their seeds and 
their capacity to fix atmospheric nitrogen due to a symbiosis 
with rhizobia, legume plants have an essential role in human 
and animal nutrition and in sustainable agriculture. Alfalfa 
(Medicago sativa) yields in Europe are strongly reduced by 
the impact of Verticillium wilt caused by V. albo-atrum. In 
alfalfa, monogenic and quantitative resistance to Verticillium 
wilt has been reported (Pennypacker, 2000; Acharya and 
Huang, 2003). Antimicrobial compounds such as phytoalex- 
ins and saponins have been suggested to be involved, as well as 
deposit of electron-dense material in xylem vessels (Acharya 
and Huang, 2003). Although tolerant cultivars are available 
(Molinero-Demilly et al, 2007), the molecular mechanisms 
underlying resistance are still not understood. This lack 
of data stresses the need for a model system to investigate 
Verticillium wilt response in legume plants. 

Because of its autotetraploid and outcrossing nature, genetic 
studies on disease resistance in alfalfa are difficult. The closely 
related wild species Medicago truncatula has been established 
during the last 1 5 years as a model for legume plants, particularly 



attractive for the study of plant-microbe interactions (Samac 
and Graham, 2007; Rose, 2008; Young and Uvardi, 2009). It is 
diploid, self-fertile, has a short generation cycle, a small genome 
size, and shows high synteny with cultivated legume crops. 
Hence, knowledge obtained by studying M. truncatula can be 
transferred to other legumes. The understanding of the molecu- 
lar mechanisms involved in symbiosis with rhizobia and mycor- 
rhiza has greatly advanced thanks to studies with this plant, and 
common and specific symbiotic pathways for these symbioses 
have been identified (Gough and Cullimore, 201 1). 

In addition to the tremendous amount of data on symbiotic 
interactions, pathogenic interactions have also been described. 
Resistance genes and loci against various pathogens and pests, 
such as Colletotrichum trifolii, Erysiphe pisi, Aphanomyces 
euteiches, Ralstonia solanacearum, and pea aphid have been 
identified in M. truncatula (Vailleau et al, 2007; Yang et al, 
2007; Ameline-Torregrosa et al, 2008; Djebali et al, 2009; 
Pilet-Nayel et al, 2009; Stewart et al, 2009). Furthermore, the 
RCT1 gene responsible for resistance against C. trifolii race 1 
in M. truncatula was able to confer resistance against several 
races of the fungus in alfalfa (Yang et al, 2008b). Yet, com- 
pared to other plant families, genetic mechanisms of disease 
resistance in legumes, particularly for root diseases, are less 
understood. They may present particular characteristics due to 
their putative interactions with symbiotic processes (Mithofer, 
2002), which are absent in non-legume plants. Such putative 
interaction with symbiosis may lead to novel recommendations 
in terms of plant breeding for root disease tolerance in legumes. 

In the present work, the model plant M. truncatula is 
used to study Verticillium wilt in legumes. The infection was 
characterized and colonization of the host plant was stud- 
ied with a green fluorescent protein (GFP)-expressing strain. 
Resistance and susceptibility levels of different M. truncatula 
lines were quantified through disease symptoms modelling. 
Loci responsible for resistance were identified in three crosses, 
and a collection of accessions was used to analyse biodiver- 
sity of the interaction. In order to address the question of 
crosstalk between disease response and symbiosis, various 
nodulation mutants in line A17 were inoculated with the fun- 
gus and symptom development was analysed. 

Materials and methods 

Plants 

M. truncatula Gaertn. seeds were obtained from plants grown in the 
greenhouse. The lines, except Jemalong-A17, and mutants are derived 
from natural populations of different geographical origins (http://wwwl . 
montpellier.inra.fr/BRC-MTRymauguio/mauguio.php). Scarification 
and germination of seeds, as well as culture in peat substrate (Jiffy pots, 
Jiffy, France), were performed as described by Vailleau et al. (2007). 
For microscopy, plants were grown in a miniaturized hydroponic culture 
system on N-containing Farhaeus medium (Farhaeus, 1957). In vitro 
plants were grown on solid Farhaeus medium as described by Vailleau 
et al. (2007). Plants were grown in a growth chamber with 16 h light 
(170 umol nT 2 s _I ) at 25 °C and 8 h dark at 23 °C. 

Fungal isolates 

V. albo-atrum V31-2 was grown on PDA medium at 24 °C in the dark. 
Spore suspensions were obtained by flooding 2-week-old cultures in 
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Petri dishes with sterile water. The spore concentration was determined 
with a Malassez haematocytometer. Wild type and GFP-expressing 
strains were maintained on PDA plates at 14 °C or stored as spore sus- 
pensions in 25% glycerol at -80 °C. 

Genetic transformation of V. albo-atrum 

V. albo-atrum V3 1 -2 was transformed with the GFPgene by Agrobacterium 
tumefaciens-mediated transformation as described by Eynck et al. (2007), 
using A. tumefaciens AGL-1 containing the binary vector pBin-GFP-hph. 
This vector carries the GFP gene under the control of the constitutive glyc- 
eraldehyde-3-phosphate dehydrogenase promoter of Aspergillus nidulans 
and a hygromycin B resistance marker (O'Connell et al, 2004). Equal 
volumes of fresh bacterial and spore suspensions were mixed and plated 
onto sterile cellophane sheets placed on co-cultivation medium containing 
200 |iM acetosyringone. After 2 days of incubation at 21 °C in the dark, 
the fungal and bacterial culture growing on the cellophane sheets was 
resuspended in 3 ml of sterile 0.9% NaCl containing 200 uM cefotaxime, 
and the suspension was plated on selective Czapek Dox medium contain- 
ing 200 |iM cefotaxime and 50 |ig.mr' hygromycin B for selection of 
transformants. Colonies of the fungus appeared after 5 days of incubation 
at 21°C and putative transformants were transferred to fresh Czapek Dox 
medium containing 50 ng.ml hygromycin B. Monosporic strains were 
obtained after sporulation of fluorescent isolates on this medium. Their 
stability was tested by four generations of culture on PDA without selec- 
tion pressure, then culture on PDA with 50 ng.ml hygromycin B and 
assessment of their fluorescence under the microscope. 

Inoculation and symptom scoring 

Inoculation experiments were performed at three sites using the same 
standard protocol: Ecolab at Auzeville, R2n at Druelle, and Barenbrug 
Tourneur at Mas Grenier. Suspensions of conidia were obtained as 
described above, and the concentration was adjusted to 10 6 spores mr 1 
for all inoculations. Ten-day-old plants (first pair of leaves stage) grown 
in peat substrate were inoculated after cutting 1 cm of the bottom of 
the Jiffy pots which were then dipped in the conidial suspension during 
30 minutes before being transferred to trays containing moist soil. For 
microscopy studies, inoculation was performed with 10-day-old plants 
grown in hydroponic conditions by cutting 1 cm of their roots and dipping 
them in conidial suspension for 30 min before transferring them back to 
the nutritive solution. For in vitro experiments, 4-day-old plants grown 
on solid Farhaeus medium were inoculated with 100 ul of conidia solu- 
tion after cutting 1 cm of the root. In some assays, 3-4-week-old plants 
grown in Jiffy pots were inoculated with conidial suspension sprayed 
on cut aerial parts as described by Molinero-Demilly et al. (2007). All 
inoculated plants and respective controls treated with water were incu- 
bated in a growth chamber at 20 °C with 1 6 h photoperiod. Symptoms 
were rated on a scale from 0 to 4, similar to that described for alfalfa 
(http://www.naaic.org/stdtests/verticil.htm, Fig. ID) and the proportion 
of dead plants (PDP) was determined at the end of the experiment. 

The area under the disease progress curve (AUDPC, Shaner and 
Finney, 1977) was computed based on the disease scores from 0 to 4, 
using the 'agricolae' package of the R system (R Core Team, 2012). 
Model for disease progress empirical curves was based on a logistic 
model. The curves were fitted to disease index data based on the over- 
all methodology described in Gilligan (1990). Model fitting was done 
using the 'nlme' package for the R system (Pinheiro et al, 2012). Full 
models had fixed components for line and when required, treatment and 
repetition effects and had random components for plants within line. 
When needed to improve model fit, autocorrelation among residuals 
was specified. Models were fitted using maximum-likelihood method, 
and analysis of variance (ANOVA) allowed testing for significance 
of effects. BLUE and BLUP of effects were obtained using the most 
parsimonious model for a given experimental design. The modelling 
allowed the description of the disease progress in terms of three func- 
tional parameters: (i) maximum disease index (MDI) at the end of the 
experiment; (ii) time to reach 50% MDI; and (iii) time to proceed from 
50% MDI to (50% MDI)/(1 + e" 1 ) (c.75% MDI). 



Microscopic observations 

Roots were rinsed briefly in water and 1-cm long fragments were 
embedded in 5% low melting point agarose. Longitudinal and radial 
sections of 80 and 45 um, respectively, were prepared on a vibratom 
(VT 1000S, Leica, Germany) and mounted on a glass slide with distilled 
water. Confocal images were acquired with a spectral confocal laser 
scanning system (SP2 SE, Leica) equipped with an upright microscope 
(DM 6000, Leica). Observations were made using xlO (HC PL Fluotar, 
N.A. 0.3) and x40 (HCX PL APO, N.A. 0.8) dry and water immersion 
objectives, respectively. The 488-nm ray line of an argon laser was used 
to detect the GFP fluorescence emission collected in the range between 
490 and 540 nm. 



Genetic maps 

Previously described low-density genetic maps of A17 x DZA315.16 
(LR4; Julier et al, 2007) and A17 x F83005.5 (LR5; Arraouadi et al., 
2012) recombinant inbred line (RIL) populations were improved by 
adding simple sequence repeat markers based on the M. truncatula 
genome sequence assembly (Young et al., 2011, http://www.medicago- 
hapmap.org/?genome) and unigene set of the Medicago Gene Index at 
DFCI (http://compbio.dfci.harvard.edu/). PCR conditions, gel electro- 
phoresis, and genotype scoring were done as previously reported (Julier 
et al., 2007). The genetic map of the F83005.5 x DZA45.5 cross (LR3) 
is described in Hamon et al. (2010). 

Quantitative trait locus analysis 

Response to V. albo-atrum was evaluated in three different RIL popula- 
tions and in different sites and/or years. For each RIL, 8-12 plants were 
inoculated as described above and scored regularly during 4 weeks. 
A total of 173 RILs of the Al 7 x F83005.5 (LR5) population in the 
F8 and/or F9 generations were evaluated in two different sites (ENSAT 
and Barenbrug-Tourneur), and AUDPC, disease curve parameters, and 
PDP were computed. Similarly, 137 RILs of the A17 x DZA315.16 
(LR4) population in the F8 and/or F9 generations were evaluated in 
three sites (ENSAT, Barenbrug-Tourneur, and R2n) and 178 RILs of the 
F83005.5 x DZA45.5 (LR3; Hamon et al., 2010) RIL populations were 
evaluated at R2n in two independent experiments and mean PDP were 
determined. 

Quantitative trait loci (QTLs) for resistance to V. albo-atrum were 
detected by Multiple QTL Mapping (MQM) (Jansen, 1993; Jansen and 
Stam, 1994) using the 'qtP package (Broman et al, 2003; Arends et al, 
2010) of the R system. AUDPC, PDP, and disease parameters were 
used as variables for QTL detection. Empirical threshold values for the 
logarithm of odd scores (LODs) were determined by computing 5000 
permutations (Churchill and Doerge, 1994). Depending on the pheno- 
type and the cross, the critical LOD score to indicate QTL significance 
ranges from 2.5 to 3.7. Heritabilities were computed using variance- 
components methods, by equating mean squares to their expectations 
and using WLS method for variance estimations (Kearsey and Pooni, 
1996). All computations were done using the R system. 

Results 

V. albo-atrum is a pathogen of M. truncatula 

To assess if the model legume M. truncatula is also a host 
plant for the alfalfa pathogen V. albo-atrum, the alfalfa iso- 
late V31-2 was inoculated on six M. truncatula genotypes 
(Jemalong A17, A20, DZA45.5, DZA315.16, F83005.5, and 
TNI .11) which are parental lines of mapped RIL populations. 

Young plants cultivated in peat substrate were root inoc- 
ulated and compatible and incompatible interactions were 
identified: pronounced wilting symptoms and death of the 
plants were observed in some lines such as F83005.5 and 
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Fig. 1. Verticillium wilt symptoms in Medicago truncatula. (A-C) Phenotypes of the two most contrasted lines of M. truncatula 
(Jemalong-A1 7 and F83005.5) after root inoculation with a conidia suspension of Verticillium albo-atrum V31-2 (10 6 spores mr 1 ) in 
peat substrate (A), hydroponic system (B), and in vitro culture (C). (D) Symptom scale for disease index scoring. Scores from 0 to 4 
correspond to characteristic stages of disease development. 



DZA315.16 whereas other genotypes like A17 and DZA45.5 
seemed to be highly resistant. 

The two most contrasting lines A 17 and F83005.5 (Fig. 1A) 
were also root inoculated in hydroponic and in vitro culture 
conditions. Again, plants developed typical disease symp- 
toms and the same differential behaviour (Fig. IB and 1C, 
respectively). 

To quantify the degree of susceptibility, symptom develop- 
ment in the six parental lines was regularly scored on a scale 
from 0 to 4 (Fig. ID). AUDPC and PDP at the end of the 
experiment were computed. Both disease parameters showed 
significant differences among the six parental lines (ANOVA 
P-values <2.10~ 16 ). The symptom curves were also modelled 
with a non-linear mixed model which allows characteriz- 
ing the disease development more precisely with three func- 
tional parameters: the MDI which measures the severity of 
the disease, the time to reach 50% MDI which describes the 
onset of disease, and the time to proceed from 50 to 75% of 
the MDI which represents the rate of disease progression 
(Gilligan, 1990) (Fig. 2A). MDI was the only parameter of 
this model which showed significant differences among the 



six lines (ANOVA P- value <0.0001). As assessed by MDI 
and AUDPC values, line F83005.5 was the most susceptible 
one, followed by DZA315.16; A20 was intermediate, and 
DZA45.5, TNI. 11, and A17 were resistant, with A17 being 
highly resistant (Fig. 2). When considering PDP, DZA45.5, 
TN1.1 1, and A17 all appeared highly resistant (Fig. 2B). Thus 
PDP seems to be less discriminating than AUDPC or MDI to 
describe disease resistance levels among genotypes. The same 
six lines were also spray inoculated on wounded aerial parts, 
as described for resistance screening in M. sativa (Molinero- 
Demilly et al, 2007), and displayed similar phenotypes as 
in the root inoculation assay (Supplementary Fig. SI, avail- 
able at JXB online). Taken together, the results showed that 
M. truncatula is a host plant for V. albo-atrum and that sources 
of resistance exist in the germplasm of this model plant. 

V. albo-atrum colonizes the host plant through xylem 
vessels 

To study colonization by the pathogen, V. albo-atrum V31-2 
was transformed, via A. tumefaciens, with the GFP reporter 
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Fig. 2. Evolution of wilting symptoms and disease evaluation parameters of six Medicago truncatula lines with different levels of 
quantitative resistance. (A) Fitted logistic curves for disease kinetics for six genotypes (F83005.5, DZA315.16, A20, DZA45.5, TN1 .1 1 , 
Jemalong-A1 7) in three independent experiments (RB1 , RB2, and RB3), following root inoculation in peat substrate with Verticillium 
albo-atrum V31-2 at 10 6 spores ml -1 . Wilting symptoms were scored for 20-30 days after inoculation, on a scale from 0 to 4. Raw 
disease curves were modelled with a non-linear mixed model. Lines are ordered by decreasing susceptibility from left to right. (B) 
Distribution of the disease evaluation parameters after root inoculation among the six lines. Lines with the same letter belong to the 
same group of means, defined by Newman-Keuls test. AUDPC, area under disease progress curve. Data were obtained in three 
independent experiments performed at INP-ENSAT and R2n sites with 8-10 plants per line per experiment. 
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gene following a protocol described for V. dahliae (Eynck 
et at, 2007). Four monosporic fluorescent strains were 
obtained. Their growth rate on PDA medium and their 
aggressiveness on the susceptible M. truncatula line F83005.5 
were similar to the wild type, and they were also unable to 
cause disease in the resistant line A17 (Supplementary Fig. 
S2). V. albo-atrum Alb2, which has high fluorescence, was 
then used for root inoculation in hydroponic cultures of the 
susceptible line F83005.5 and for microscopic observations. 
Roots were harvested at various times after inoculation, and 
sections were observed by confocal laser scanning micros- 
copy. Conidia were rapidly absorbed into the xylem vessels 
by the transpiration stream as seen at 2 hours post inocu- 
lation (Fig. 3A). It was observed that their movement was 
stopped by the plaques which could then be crossed by ger- 
minating hyphae (Fig. 3B). Germination occurred as soon as 
after 3 hours. Hyphae then grew inside the vessels along the 
root axis but were also able to penetrate into adjacent xylem 
vessels (Fig. 3C). Sporulation was observed in xylem vessels 
after several days (data not shown). After 12 days, hyphae 
were detected in the hypocotyl, but the fungal growth was still 
restricted to the central cylinder. At later stages, when plants 
were severely diseased (i.e. symptom score >3.5 after 21 days), 
a massive production of mycelium was observed and hyphae 
colonized stele parenchyma, penetrated through the endoder- 
mis, and colonized the root cortex (Fig. 3F and 3G). 

Taken together, the colonization pattern of V. albo-atrum 
in M. truncatula was typical for other described vascular wilt 
pathogens (Agrios, 2005; Fradin and Thomma, 2006) and con- 
firmed that V. albo-atrum is a true pathogen of M. truncatula. 

Evidence for diverse quantitative genetic mechanisms 
controlling resistance to V. albo-atrum in M. truncatula 

Resistance parameters segregate in different RIL 
populations 

To search for loci controlling the resistance to V. albo-atrum, 
RILs derived from crosses involving the most contrasting 
parental lines were analysed. Disease symptoms were 
screened within three RIL populations, namely 178 RILs of 
the DZA45.5 x F83005.5 (LR3) population, 137 RILs of the 
A17 x DZA315.16 (LR4) population, and 173 RILs of the 
A17 x F83005.5 (LR5) population. 

The PDP at the end of the experiment was evaluated in 
different sites for all three populations (21-30 days post inoc- 
ulation, depending on the site). For LR4 and LR5 popula- 
tions, foliar wilting symptoms were also scored regularly and 
AUDPC and disease functional parameters were computed 
(Supplementary Fig. S3). For these two crosses, three inde- 
pendent experiments were performed in at least two differ- 
ent sites. LR3 population was analysed in two independent 
blocks at one site and mean PDP values were determined. 
ANOVA showed a significant effect of the trial (block or site) 
for most crosses and variables (Table 1). Yet, the high positive 
Spearman's correlations obtained on genotype ranks for dif- 
ferent repetitions showed that relative susceptibility levels of 
the RILs were reproducible over the different environments. 



Among the four parameters describing the disease kinetics 
in LR4 and LR5 RILs, only MDI and AUDPC values were 
significantly different between genotypes (Table 1). ANOVA 
also showed significant genetic variation (P < 2.10~ 16 ) for the 
PDP at the end of the experiment in the three studied RIL 
populations. As expected, highly significant positive correla- 
tions were observed between these three parameters (Pearson's 
correlation coefficient R > 0.95 in each population). 

For each of these traits, some RILs showed lower or higher 
mean values compared to parental means, suggesting trans- 
gressions in phenotypes (Supplementary Fig. S4). Broad 
sense heritabilities, that describe the part of total phenotypic 
variation due to genetic variation, range from 31 to 33% for 
AUDPC and MDI values and can reach up to 65% for the 
PDP in LR5 population (Table 2). Taken together, the exist- 
ence of transgressive lines and moderate to high heritabilities 
for all phenotypes confirm that the populations and traits are 
suited to QTL mapping through integration of RILs geno- 
typing and phenotyping data. 

Two distinct combinations ofQTLs on chromosome 7 and 
chromosomes 2 and 6 govern resistance to Verticillium. 

To facilitate precise QTL localization, marker densities were 
increased for LR4 (Julier et al, 2007) and LR5 (Arraouadi 
et al, 2012) genetic maps. Additional simple sequence repeats 
were thus genotyped in these two RIL populations. The LR4 
map now comprises 370 loci, evenly distributed over each link- 
age group (Supplementary Fig. S5A). It spans a total of 910 
cM with an average distance between markers of 3 cM. The 
LR5 map comprises 1 58 loci and spans a total of 780 cM with 
an average distance between markers of 5 cM (Supplementary 
Fig. S5B). The linkage groups of the genetic maps presented 
here are longer than in previously published maps. This shows 
that detection of crossing-over events was more accurate by 
increasing the number of mapped simple sequence repeats. 

The RIL populations LR4 and LR5, both derived from 
crosses involving the resistant Jemalong-A17 line as a par- 
ent, were first analysed for response to V. albo-atrum. A sin- 
gle QTL, located on linkage group 7, was detected for both 
populations and with each of the three disease parameters 
(Table 2). The 1-LOD support intervals of the QTL for 
each trait ranged from 7 to 15 cM and were overlapping. 
This suggests that, for both populations, the same genomic 
region between markers MTE73 and MTE126, is involved 
in the genetic control of Verticillium wilt (Fig. 4). This single 
major QTL, named Mt Val, was detected with data obtained 
in different sites for both populations (Supplementary Fig. 
S6). The Mt Val QTL showed high LOD scores (up to 26 
for AUDPC in LR5) and accounted for large effects, as it 
explained from 20.6% of the phenotype for the PDP in LR4 
to 42.5% of phenotypic variation for AUDPC parameter in 
LR5. The occurrence of one single QTL suggests also that 
the two susceptible lines F83005.5 and DZA315.16 do not 
contain additional resistance loci with strong effects. 

In order to study resistance in a parent different from 
Jemalong-A17, the cross F83005.5 (susceptible) x DZA45.5 
(resistant) (LR3) was used. QTL analysis of PDP in this 
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Fig. 3. Root colonization of Medicago truncatula F83005.5 (susceptible line) by a GFP-expressing strain of Verticillium albo-atrum, 
observed by confocal laser scanning microscopy. (A) At 2 hours post inoculation (HPI), Conidia are seen in xylem vessels and 
accumulate at the end of a tracheid (arrow). (B) Conidia penetrating with its germination hypha through the plaque at the end of a xylem 
vessel (enlargement of A). (C) At 24 HPI, germinating conidia and hyphae. (D) At 7 days post inoculation (DPI), fungal hyphae in the 
central cylinder; no hyphae detected in the cortex. (E) At 7 DPI, proliferation of mycelium in xylem vessels and penetration into adjacent 
vessels (asterisk). (F and G) At 21 DPI, proliferation of mycelium in vascular (F) and cortical (G) tissues. Bars = 50 jim, and 10 jim in 
B. Longitudinal sections in A, B, C, E, G; cross-sections in D and F. cor = cortex; cw = cell wall; hy = hypha; sp = conidia; xy = xylem 
elements. 



mapping population evidences one major locus on chromo- 
some 2 and another minor locus on chromosome 6 (Table 2). 
These QTLs, called MtVa2 and MtVa3, had respective LOD 
scores of 6.0 and 3.6, and together accounted for 22% of the 



resistance phenotype and covered confidence intervals of 4 
and 21 cM. 

At each of the three loci detected by analysis of the differ- 
ent RIL populations, the allele from the resistant parent (A17 
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Table 1 . Variance analysis for several disease evaluation parameters in three RIL populations: proportion of dead plants at the end of the 
experiment, disease functional parameters, and AUDPC Significance levels are indicated for each factor. For LR3 population, only proportion 
of dead plants at the end of the experiment is available. [R] and [S] indicate resistant and susceptible parental lines, respectively. AUDPC, 
area under the disease progress curve; NS, not significant at a = 1 %. Percentage of dead plants was tested by generalized linear model. 



Cross 



Source of 
variation 



df 



Percentage of 
dead plants 



Disease curve modelling 



Maximum 
disease index 



Time to 50% 
disease (days) 



Rate of disease 
progression (days) 



AUDPC 



LR3: F83005.5 [S] x DZA45.5 [R] Block 

Genotype 



1 NS 

212 <2.1CT 16 



LR4:A17[R]xDZA315.16[S] Site 2 2.10" 

Genotype 158 <2.1CT 



<2.1CT 16 
4.10- 10 



<2.icr 16 

NS 



NS 
NS 



<2.icr 16 
4.1 cr 9 



LR5: A1 7 [R] x F83005.5 [S] Site 1 NS 

Genotype 159 <2.1CT 



<2.icr 

1.10- 



<2.10- 16 
NS 



NS 
NS 



3.10" 
1.10-' 



Table 2. Characteristics for the quantitative trait loci of resistance to Verticillium albo-atrum detected in LR3, LR4, and LR5 recombinant 
inbred line populations [R] and [S] indicate resistant and susceptible parental lines, respectively, h 2 , broad sense heritability; (H), narrow 
sense heritability; LOD, logarithm of odd score; QTL, quantitative trait locus. 



Cross 


Trait 


Trait 


Chromosome 


Nearest 


Peak 


Peak 


1-LOD 


QTL 


Percentage 


Percentage 






heritability 




marker 


LOD 


(cM) 


support 


effect 


of phenotypic of additive 






h 2 (H) 






score 




interval (cM) 




variance 


variance 


LR3: F83005.5 [S] x 


Percentage of 


0.32 (0.19) 


2 


EM4.250 


6.0 


196 


193-197 


-0.095 


22.0 


36 


DZA45.5 [R] 


dead plants 




6 


mtic343 


3.6 


142 


133-154 


-0.067 




18 




Maximum 


0.33 (0.20) 


7 


mtic432 


11.1 


60 


55-70 


0.41 


32.2 


92 




disease index 




















LR4: A1 7 [R] x 


AUDPC 


0.31 (0.18) 


7 


mtic273 


12.3 


58 


55-70 


0.0022 


32.9 


83 


DZA315.16 [S] 


Percentage of 


0.46 (0.29) 


7 


mtic82 


5.7 


76 


65-78 


0.0907 


20.6 


32 




dead plants 






















Maximum 


0.32 (0.19) 


7 


MTE85 


22.6 


58 


54-61 


0.5093 


39.4 


99 




disease index 




















LR5: A1 7 [R] x 


AUDPC 


0.32 (0.19) 


7 


MTE85 


26.1 


58 


54-61 


0.0028 


42.5 


99 


F83005.5 [S] 


Percentage of 


0.65 (0.48) 


7 


MTE85 


14.9 


65 


56-68 


0.1973 


29.3 


48 




dead plants 





















or DZA45.5 depending on the concerned RIL populations) 
was associated with lower susceptibility level computed as 
AUDPC, MDI, or PDP (Table 2). 

V. albo-atrum resistance QTLs do not co-localize with 
identified Ve homologues 

To check if Ve gene homologues are localized in the region 
underlying the QTLs, a blastp analysis on the M. truncatula 
genome was performed using the peptide sequence of tomato 
VEX (GeneBank accession AAK58681.2). Eight homologous 
sequences were detected which were further analysed by mul- 
tiple alignment and establishment of a phylogenetic tree using 
clustal w version 2 (http://www.ebi.ac.uk/Tools/msa/clustalw2/) 
(Supplementary Fig. S7A). Among them, seven are localized on 
chromosome 4 and one on chromosome 5. Hence these homolo- 
gous sequences do not co-localize with the major QTLs of par- 
tial resistance to V. albo-atrum in this study's pathosystem. 

The four most homologous candidates (Medtr4g017280.1, 
Medtr4g017350.1 Medtr4g017370.1, and Medtr5g046350.1) 



were compared with the complete open reading frame amino- 
acid sequences of Ve genes described to be involved in resist- 
ance to Verticillium wilt, such as Vel and Ve2 from tomato, 
StVe from Solanum torvum, GbVe from cotton, as well as 
Ms Vel and MWel from Mentha spicata and Mentha latifo- 
lia respectively (Supplementary Fig. S7B). The comparison 
shows a high diversity in the peptide sequence of the vali- 
dated Ve proteins with identities ranging from 42 to 96%. The 
four M. truncatula homologues share from 69 to 91% iden- 
tity. The highest homology between the candidates and vali- 
dated Ve proteins was 46% identity for three Mt genes with 
Vel. The four candidates have similar protein structures as Ve, 
such as LRRs, signal peptides, and transmembrane domains 
(Supplementary Fig. S7C). 

M. truncatula presents a wide biodiversity in its 
response to V. albo-atrum 

The QTL analysis revealed that different resistance sources 
exist within the species. To gain a broader view of this finding, 
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Fig. 4. Position of quantitative trait loci (QTL) for Verticillium albo-atrum resistance in Medicago truncatula assessed through maximum 
fitted disease index and AUDPC in LR4 and LR5 populations and proportion of dead plants in LR3, LR4, and LR5 populations. Bars at 
the left of considered linkage groups indicate 1 logarithm of odd score support for QTL interval. 



the natural biodiversity in response to V. albo-atrum infection 
was analysed by phenotypic evaluation of a M. truncatula 
core collection of 32 lines (CC32), representing the genetic 
variability within this Mediterranean natural weed (Ronfort 
et al, 2006). ANOVA indicated statistically significant differ- 
ences in terms of disease onset and severity in the responses 
of these lines (Fig. 5A, Table 3, and Supplementary Table SI). 
The response to V. albo-atrum of the CC32 lines appeared 
to be independent of the M. truncatula subpopulation struc- 
ture defined with neutral markers by Ronfort et al. (2006). 
This result (Fig. 5B), was confirmed through ANOVA, which 
showed that P-values for the subpopulations effect were 
non-significant for MDI and for disease onset (P - 0.26 and 
P - 0.12, respectively). 

Phenotypic screening of the CC32 also showed that the 
speed of disease onset was independent of the symptom 
severity at the end of the experiment (Fig. 5B). Thus, classes 
of typical behaviour could be identified. Among highly sus- 
ceptible lines, a group represented by DZA315.16, presented 
fast appearance of symptoms whereas another group, typi- 
fied by L321, showed a slower evolution of the disease. The 
same variability of response kinetics was observed in the most 
resistant lines, with some lines such as A17 that showed weak 
symptoms even early after inoculation and some genotypes 



such as L369 for which weak symptoms developed late. These 
results suggest the existence of diverse genetic systems con- 
trolling the resistance to Verticillium within M. truncatula 
species. 

This evaluation is leading to the identification of even 
more contrasted lines than the parents used in the present 
study, notably more susceptible ones which exhibit fast and 
intensely evolving wilting symptoms. In the future, it could be 
interesting to study the progenies of new crosses involving the 
most contrasted lines; such populations are currently under 
development. 

Response to Verticillium interacts with nodulation 
pathways 

As a legume plant, M. truncatula has the ability to inter- 
act with both root symbiotic or pathogenic microbes. As a 
preliminary survey to identify putative crosstalks in signal- 
ling pathways involved in plant defence and nodulation, ten 
mutants of M. truncatula known to be affected in distinct 
steps of the rhizobia symbiosis were inoculated with V. albo- 
atrum (for review, see Kouchi et al, 2010). NFP and DMI2 
membrane receptors perceive the Nod factors released by the 
soil bacterium Sinorhizobium meliloti and initiate the signal 
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Fig. 5. Diversity in the response of Medicago truncatula to Verticillium albo-atrum within the Core Collection 32. (A) Fitted logistic 
curves for disease kinetics; M. truncatula lines are ordered from top left to bottom right by increasing maximum disease index, from the 
resistant line Jemalong-A1 7 (L738) to the most susceptible line L321 ; similar results were obtained in two independent experiments. (B) 
Ouantitative resistance levels as evaluated through maximum fitted disease index and time to reach 50% of maximum disease index, 
adjusted over two biological repeats. Each of the four putative M. truncatula sub-populations (Ronfort ef al. 2006) is identified by a 
distinct symbol (1 =9 2=M, 3=*, 4= A). DZA45.6 is a sister line of DZA45.5. 



cascade involving proteins such as DMI1, DMI3, and NSP1 
to elicit the gene expression programme leading to infection 
and nodulation. HCL, LIN, and RPG play roles in the rhizo- 
bia infection and SUNN and SICKLE participate to the 
autoregulation of the symbiotic process. These mutants are 
in the genetic background of the resistant line Jemalong-A17. 

The hypernodulation mutants sickle and sunn showed 
altered but opposed responses to V. albo-atrum (Table 4 and 
Supplementary Fig. S8). As assessed by all three parameters 
of the disease curve, sunn was significantly more susceptible 
than the wild type, with increased wilting symptoms and 
faster disease development. In contrast, sickle seemed to be 
more resistant than the wild type as evidenced by the clearly 
slower disease progression and the subsignificant reduction 
in symptom severity. The nfp mutant appeared to be more 
susceptible to V. albo-atrum than the wild type, as revealed 
by the slight increase in MDI. Besides these three mutants 
that showed highly altered disease phenotypes, three other 
ones, namely dmil, lin, and nspl, also showed some increase 
in their resistance level as revealed by the significant reduc- 
tion in the disease progression rate. 



Discussion 

In order to study Verticillium wilt in legumes, a new patho- 
system involving the model legume M. truncatula and an 
alfalfa isolate of V. albo-atrum was established in this study. 
Susceptible and resistant lines were identified with differ- 
ent inoculation methods that reproducibly established lines 



A17 and DZA45.5 as resistant and lines F83005.5 and 
DZA315.16 as susceptible. The colonization of susceptible 
roots as observed with a GFP-expressing V. albo-atrum strain 
was limited to the vascular system until the final stage of dis- 
ease when the fungus invaded the whole plant tissue. This is in 
agreement with colonization patterns reported in other host 
plants (Fradin and Thomma, 2006) and confirms M. trunca- 
tula as a true host plant for V. albo-atrum. The tools available 
for this model plant allowed the study of the genetic basis of 
resistance to V. albo-atrum. 

Analysis of three RIL populations led to the identification 
of three QTLs for resistance in two resistant parental lines. 
These loci are located on linkage group 7 for A 17 and on 
linkage groups 2 and 6 for DZA45.5, respectively. The analy- 
sis revealed thus that only a small number of major QTLs 
are involved in resistance to V. albo-atrum in M. truncatula, 
suggesting a simple genetic control. Model fitting of disease 
progression curves in RIL populations allowed this study to 
show that these resistance traits are more characterized by 
lower symptoms severity at the end of the experiment, rather 
than delay in disease onset or lower disease progression rate. 
AUDPC combines all the three functional parameters and 
may be less powerful in discriminating between resistance 
QTLs for disease progression and resistance QTLs for symp- 
toms severity. These two types of resistance traits probably 
involve distinct genetic mechanisms. 

Experiments in three different sites with two RIL popu- 
lations using A17 as the resistant parent led to the identi- 
fication of the same major QTL on linkage group 7. Thus, 
while environmental variations significantly affect the RILs' 
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Table 3. Disease functional parameters in the core collection 32 of Medicago truncatula. Lines are ordered in increasing maximum 
disease index parameter (i.e. from most resistant to most susceptible. Plants were grown in peat substrate and inoculated by a 
conidial suspension of Verticillium albo-atrum (10 6 spores mr 1 ). Disease symptom scoring was done in two independent experiments. 
Populations are noted according to Ronfort etal. (2006). CC32, core collection 32. 



CC32 line code 


Usual name 


Maximum disease 


Time to 50% disease 


Geographical 


Population 






index 


(days) 


origin 




L245 


SA014161 


1.26 


15.17 


Jordan 


1 


L738 


Jemalong A1 7 


1.29 


14.51 


Australia 




L736 


DZA045.6 


1.58 


15.14 


Alqeria 


3 


L369 


PRT1 80-A 


1.60 


18.64 


Portuqal 


2 


L544 


ESP105-L 


1.78 


17.93 


Spain 


2 


L648 


SALSES 42B 


1.86 


17.39 


France 


2 


L542 


DZA233.4 


1.92 


14.55 


Alqeria 

.J 


3 


L1 44 


SA014163 


1.93 


16.56 


Jordan 


3 


L174 


SA028064 


1.97 


14.90 


Cyprus 


1 


L545 


ESP 158-A 


2.03 


13.87 


Spain 


2 


L679 


F6601 7 


2.04 


16.26 


France 


2 


L554 


F20089-B 


2.09 


17.45 


France 


3 


L213 


SA027882 


2.17 


14.69 


Morocco 


4 


L550 


F1 101 3-3 


2.25 


19.16 


France 


4 


L310 


SA009944 


2.30 


14.98 


Tunisia 


3 


L290 


SA009119 


2.32 


15.63 


Turkey 


3 


L337 


GRC 042-1 


2.38 


16.12 


Greece 


1 


L530 


F83005-5 


2.51 


15.05 


France 


3 


L368 


DZA012-J 


2.63 


16.06 


Algeria 


4 


L543 


DZA 327-7 


2.64 


16.07 


Algeria 


2 


L239 


SA026063 


2.69 


15.61 


Morocco 


2 


L154 


SA024714 


2.70 


14.71 


Italy 


1 


L555 


GRC 020-B 


2.74 


15.87 


Greece 


1 


L734 


DZA315.16 


2.83 


13.70 


Algeria 


3 


L049 


SA009707 


2.87 


15.40 


Tunisia 


4 


L557 


GRC 064-B 


2.88 


15.60 


Greece 


3 


L163 


SA022322 


2.91 


15.62 


Syria 


1 


L651 


SALSES 71 B 


3.02 


17.29 


France 


4 


L549 


F11005-E 


3.03 


17.12 


France 


4 


L552 


F20047-A 


3.18 


15.08 


France 


1 


L198 


SA009048 


3.30 


13.86 


Libya 


3 


L263 


SA003116 


3.38 


15.14 


Israel 


3 


L321 


SA003780 


3.57 


17.08 


Italy 


1 



responses to Verticillium (Table 1), they do not modify 
the genetic control for resistance. Such a stability in QTL 
effects towards environments and resistance parameters has 
already been reported in several pathosystems (Buerstmayr 
et al, 2009; Hamon et al, 2011; Lu et al, 2012). A survey 
of the biodiversity of the response to V. albo-atrum infection 
revealed a large variety of response, in terms of disease onset 
and severity, suggesting that different genetic systems may 
exist within the species. Response to Verticillium wilt does not 
appear to be correlated with population structure based on 
neutral markers (Ronfort et al, 2006) which represents the 
ancient history of the species. This suggests that resistance 
to Verticillium is being selected within environments. In their 
native environment, M. truncatula populations may encoun- 
ter different Verticillium strains, against which distinct resist- 
ance mechanisms may have evolved. Depending on the QTL, 
the resistance towards V. albo-atrum V31-2 may or may not 
be affected by environmental conditions as exemplified by the 
range of additive variance explained by this QTL (36% for 



the LR4 and 99% for LR5 RILs populations). Thus, whether 
or not the resistance mechanisms evolved in the natural pop- 
ulations are influenced by environmental variations remains 
to be studied. A more in-depth assessment of this hypoth- 
esis may be promoted by mapping QTLs in additional RIL 
populations derived from crosses between lines with strong 
contrasts in the response to Verticillium. Additionally, puta- 
tive correlation of the resistance with geographical origin of 
the accessions may be detected using a wider collection of 
genotypes and comparing with eco-geographical covariates. 

Genetic studies of resistance or tolerance against 
Verticillium wilt have been reported for several plant spe- 
cies (Schaible et al, 1951; Diwan et al, 1999; Bae et al, 
2008; Yang et al, 2008a). The best-characterized example is 
monogenic resistance against V. albo-atrum and V. dahliae in 
tomato, where the resistance gene Vel is effective against race 
1 of both species (Diwan et al, 1999; Fradin et al, 2009). 
Homologues of the Ve gene have been reported in other spe- 
cies as well (Fradin and Thomma, 2006). In M. truncatula, 
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Table 4. Comparison of disease functional parameters of Medicago truncatula nodulation mutants Mutants are in the Jemalong-A1 7 
(resistant) genetic background. Plants were grown in peat substrate and inoculated by a conidial suspension of V. albo-atrum (10 6 
spores mr 1 ). Disease symptom scoring was done in three independent experiments. P-values refer to test for difference between the 
parameter of the mutant line compared to that of the wild reference Jemalong-A1 7. Bold indicates significant (P < 0.05) or subsignificant 
(P < 0.10) differences compared to the reference line Jemalong-A17. [R] and [S] indicate resistant and susceptible lines, respectively. 



Line 


Maximum disease index 


Time to 50% disease (days) 


Rate of disease progression (days) 


Value 


P 


Value 


P 


Value 


P 


Mutants 














nfp 


1.3 


0.051 


24.6 


0.210 


2.7 


0.293 


dmi2 


I . I 


U.Zoo 


i y. i 


n -i a a 

U. 1 44 




u.o I u 


dmil 


1.0 


0.387 


19.7 


0.311 


3.8 


0.001 


dmi3 


1.1 


0.176 


21.2 


0.754 


1.9 


0.459 


nspl 


0.9 


0.640 


24.2 


0.244 


3.4 


0.021 


hcl 


0.6 


0.465 


21.3 


0.787 


2.4 


0.505 


Vn 


0.9 


0.680 


22.1 


0.865 


2.9 


0.064 


rpg 


1.0 


0.293 


24.2 


0.264 


1.6 


0.219 


ski 


0.4 


0.092 


24.8 


0.210 


3.2 


0.013 


sunn 


1.3 


0.042 


18.3 


0.065 


0.8 


0.000 


References 














A1 7 [R] 


0.8 




21.8 




2.2 




F8300.5 [S] 


2.8 


0.000 


16.7 


0.004 


2.0 


0.636 



four Ve homologues, exhibiting 41-46 % identity at the pro- 
tein level with Ve genes from other plant species were detected. 
As described for tomato Vel and Vel which occur in tandem 
(Fradin et at, 2009), most of the putative M. truncatula Ve 
homologues were localized near each other. Although these 
genes are on chromosomes 4 and 5 and do not co-localize 
with the resistance QTLs identified in the present study, it can- 
not be excluded that another homologue exists in the unse- 
quenced regions underlying part of the QTLs. Although 94% 
of the gene space of M. truncatula genome is now sequenced 
(Young et at, 2011), information is still lacking for some 
regions. Taking into account the large confidence interval of 
the MtVal QTL covering at least 5 Mbp and the existence 
of sequence gaps, identification of putative candidate genes 
predicted in this region was not possible. 

In Arabidopsis, overlaps between plant development 
and susceptibility to Verticillium wilt have been described. 
Notably, early flowering was reported to be correlated to sus- 
ceptibility to the Brassicaceae-infectmg Verticillium longispo- 
rum (Veronese et ai, 2003). QTLs involved in colonization by 
V. longisporum were reported to be localized in the vicinity 
of loci for the regulation of flowering such as constans and 
flowering locus (Haffner et at, 2010). A Constans-like gene 
on chromosome 7 has been proposed recently to control flow- 
ering date in M. truncatula (Pierre et at, 2011). This gene 
localizes 25 cM south of the QTL involved in resistance to 
V. albo-atrum, indicating that M. truncatula resistance does 
not seem to interfere with regulation of flowering, under the 
experimental conditions. 

M. truncatula is now well established as a model for the 
study of plant-pathogen interactions, involving bacteria, 
fungi, oomycetes, nematodes, and aphids (Torregrosa et at, 
2004; Vailleau et at, 2007; Ameline-Torregrosa et at, 2008; 
Dhandaydham et at, 2008; Djebali et at, 2009; Stewart 



et at, 2009; Uppalapati et at, 2009; Ramirez-Suero et at, 
2010). The lines A17 and F83005.5 are often compared for 
pathogenic interactions and used for comparative and/or 
genetic approaches to study the interaction with C. trifolii 
(Ameline-Torregrosa et at, 2008), A. euteiches (Djebali et at, 
2009), R. solanacearum (Vailleau et at, 2007), and Fusarium 
oxysporum (Ramirez-Suero et at, 2010). Genetic studies 
using RIL populations derived from a cross between lines 
A17 and F83005.5 revealed different quantitative mecha- 
nisms of resistance to R. solanacearum, A. euteiches, and 
C. trifolii, which always rely on a small number of loci. Two 
major QTLs against A. euteiches were reported in overlap- 
ping genomic regions of LG3 in two different RIL popula- 
tions, LR5 (Djebali et at, 2009) and LR3 (Pilet-Nayel et at, 
2009). These populations were also used by this study which, 
in contrast, led to reveal that different QTLs control resist- 
ance in the parents Jemalong-A17 and DZA45.5. So far, only 
one disease resistance gene has been cloned in M. truncatula. 
RCT1 controls resistance against C. trifolii in line A17 and 
codes for a protein of the TIR-NBS-LRR type. Introduction 
of RCT1 from M. truncatula into alfalfa conferred broad 
resistance against the anthracnose fungus (Yang et at, 
2008b), confirming the utility of model plants for knowledge 
transfer to cultivated crops. 

Compared to the model plant Arabidopsis thaliana, 
M. truncatula allows studying crosstalk between pathways 
involved in symbiotic and pathogenic interactions. Such 
crosstalk has been demonstrated by Penmetsa et at (2008) 
by using sickle, one of the numerous nodulation mutants in 
M. truncatula. This mutant was first identified by its hyper- 
nodulating phenotype and the Sickle gene was shown to be 
an orthologue of the A. thaliana Ein2 gene involved in eth- 
ylene signalling. It was also found more susceptible in its 
interactions with the root pathogens Rhizoctonia solani and 
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Phytophthora medicaginis. Four out of the ten tested nodu- 
lation mutants, ski, dmil, lin, and nspl, showed putatively 
enhanced resistance to V. albo-atrum as shown by reduced 
disease progression. However, the high resistance of the wild- 
type A17 makes enhanced resistance levels difficult to detect. 
In order to confirm these observations, a screening with 
V. albo-atrum strains for which line A17 is susceptible will be 
performed. 

In this study, sickle was slightly more resistant than the wild 
type, indicating that the ethylene signalling pathway could be 
necessary for Verticillium wilt development. The Sickle gene 
(Mtr7gl01410, corresponding to TC176805; Penmetsa et at, 
2008) is located near locus mtic337, approximately 50 cM 
south of the Mtl^7 QTL peak on chromosome 7. Data from 
the literature indicate that ethylene signalling is involved in 
the response to Verticillium wilt, but its role seems to be com- 
plex since different mutants in the ethylene pathway have dif- 
ferent disease phenotypes (Veronese et ah, 2003; Johansson 
et al, 2006; Pantelides et at, 2010). 

The second hypernodulation mutant, sunn, was also shown 
to be modified in response to V. albo-atrum in this study, but 
the phenotype was opposite to that of ski (i.e. the mutant was 
less resistant than the wild type). The sunn mutant is altered 
in a CLAVATAl-like LRR receptor-like kinase (RLK) 
(Schnabel et al, 2005). The Sunn gene (Medtr4g070970, 
which corresponds to accession AY769943) is located on link- 
age group 4. Besides hypernodulation, the mutant also shows 
increased mycorrhization (Sagan et al, 1995). Although the 
role of LRR-RLK in plant disease resistance is well estab- 
lished now, this is the first time that the sunn mutant has been 
shown to be modified in response to a pathogen. 

The nfp mutant, which is defective in a putative Nod 
factor receptor, is also affected in its resistance to V. albo- 
atrum. Since the backbone of this factor is a chitin oli- 
gomer, a modified phenotype might suggest involvement 
of PAMP perception such as chitin-derived compounds of 
fungal walls. 

As far as is known, this study is the first where a collec- 
tion of nodulation mutants was assessed for its response to a 
pathogen. Several mutants, affected in different steps of the 
nodulation process, show altered response to Verticillium wilt. 
The results emphasize the existence of crosstalks between sig- 
nalling pathways in symbiosis and disease, and the need for 
detailed studies of Verticillium wilts and other root pathogens 
in legume plants. Such crosstalks with symbiosis should be 
considered in breeding for legume tolerance for root diseases, 
in order to avoid putative trade off between nodulation and 
response to root diseases. It also raises questions concerning 
hopes to transfer rhizobial symbioses to non-legumes, as it 
may endanger the tolerance mechanisms to root disease in 
these crops. 

Future work aims at cloning the MtVal resistance gene 
by fine mapping the major QTL. However, remaining gaps 
in the sequences under this QTL render this objective a 
long-term goal relying on additional sequencing efforts in 
M. truncatula. Transcriptomic studies have now been initi- 
ated to describe the molecular response of the resistant A17 
and the susceptible F83005.5 lines. The gene expression 



profiles might also reveal candidate resistance genes if regu- 
lated at the transcriptional level. As more M. truncatula gen- 
otypes and Verticillium sp. strains are evaluated, there will 
be a focus on the discovery of new resistance loci through 
an enlargement of the analysis of biodiversity and subse- 
quent exhaustive genome-wide genetics association studies. 
This may give insights towards adaptation to root diseases in 
Mediterranean legumes. 

Supplementary material 

Supplementary data are available at JXB online. 

Supplementary Fig. SI. Comparisons of the quantitative 
resistance levels of six Medicago truncatula lines depending 
on inoculation method with Verticillium albo-atrum. 

Supplementary Fig. S2. Analysis of the growth and 
aggressiveness of four monosporic isolates of Verticillium 
albo-atrum expressing the green fluorescent protein gene in 
comparison to the wild-type V31-2. 

Supplementary Fig. S3. Example of model fitting for 
wilting symptoms in 113 LR4 recombinant inbred lines at 
site R2N (A) and 158 LR5 recombinant inbred lines at site 
ENSAT (B). 

Supplementary Fig. S4. Distribution of disease evaluation 
parameters in LR3, LR4, and LR5 recombinant inbred lines 
upon infection with Verticillium albo-atrum V31-2. 

Supplementary Fig. S5. Dense genetic maps for A17 x 
DZA315.16 (LR4, A) and A17 x F83005.5 (LR5, B) recom- 
binant inbred line populations of Medicago truncatula. 

Supplementary Fig. S6. Whole-genome scans for quanti- 
tative trait loci detection for Verticillium resistance assessed 
through maximum disease index and area under disease 
progress curve parameters in LR4 and LR5 populations for 
each site. 

Supplementary Fig. S7. Identification of four Medicago 
truncatula paralogues putatively orthologous to Verticillium 
resistance genes. 

Supplementary Fig. S8. Response to Verticillium albo-atrum 
inoculation of Medicago truncatula nodulation mutants. 

Supplementary Table SI. Variance analysis for func- 
tional disease parameters in the core collection 32, including 
Jemalong-A17. 
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